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Abstract 

The propagation of a nonlinear stationary thermomagnetic wave in superconducting media is qual- 
itatively analyzed in view of the effect of an external magnetic field. The velocity and the thickness 
of the wave front are estimated. It is demonstrated that the inclusion of an external magnetic field 
affects the thermomagnetic wave profile only slightly. 



The evolution of Joule heat during the dissipative motion of a magnetic flux gives rise to 
pronounced nonlinear effects showing up as various wave processes. These processes define 
the dynamics of thermal and electromagnetic perturbations in superconducting media [1-3] . 
The stationary propagation of a nonlinear shock thermomagnetic wave in the resistive state 
of a superconductor is an example. As is shown [1], the dispersive, nonlinear, and dissipa- 
tive processes in superconducting media result in the formation of stable thermoelectric, E, 
and/or thermomagnetic, H, waves depending on the surface conditions. It was assumed that 
the critical current density j c is independent of the local value of the external magnetic field 
H (Bean's model of critical state [4]). However, in sufficiently strong fields, this dependence 
cannot be disregarded (see, e.g., [5]) and its effect on the nonlinear thermomagnetic wave 
should be taken into account. In this paper, we studied the structure of a shock thermo- 
magnetic wave for an arbitrary dependence of the critical current density j c on the external 
magnetic field H. Expressions for the velocity and the thickness of the wave front were 
obtained. The propagation of a thermomagnetic wave in a superconducting medium along 
the x axis with a constant velocity v is described by the nonlinear heat conduction equation 
(with the self-simulated variables £(x,t) = x — vt) [1] 



— vv 



dT 
It 



d_ 



dT 



K- 



+ JE, 



Maxwell equations 



dE 



4irv 

—J, 

cr 



(1) 



(2) 



E = -H. 

c 



(3) 



and the related equation of critical state 
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j=j(T,H)+j r (E), 
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Here, v and k are the specific heat capacity and the thermal, conductivity and j r is the 
resistive current density. In the weak heating approximation (T — To) << (T c — To), the 
resistive current density j r in the range of viscous flow (E > Ef); where Ef is the effective 
field [5]) is a linear function of the vortex electric field; i.e., j r ~ crjE; where Of is the 
effective conductivity and T and T c are the equilibrium and the critical superconductor 
temperatures. For weak fields (E < Ef), the dependence j r (E) is substantially nonlinear 
and may be attributed to the thermally activated motion of the magnetic flux (flux creep 
[5]). In the range between two critical magnetic fields H Cl « H « H C2 the dependence 
of the critical current j c (T, H) on T and H can be approximated as the product of two 
functions of single variable: 



j c = j c (T,H)=j e (T)$(H). 
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Here, j c (T) = j c (To) — a(T — T ); j c (T ) = j is the equilibrium current density, a charac- 
terizes the thermal reduction of Abrikosov vortex pinning on crystal defects, and $ will be 
determined below for different values of the external field H. The thermal and electrody- 
namic boundary conditions for Eqs. (l)-(5) are given by 
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where E e is the constant external electric field. 

Jointly solving Eqs. (l)-(5) with boundary conditions (6) yields the following equation for 
the E wave: 
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is the magnetic 
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penetration depth, ad is the differential conductivity, and H e is the external magnetic field. 

Following [6], let us represent the equation for stationary points in the phase plane 
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In not-too-strong fields, when H « H C2 the $ ^— E^j function can be chosen in terms of 
the Kim- Anderson model [7]. In our situation, Eq. (10) has two equilibrium states: a stable 
node, Eq = 0, and a saddle, E = E\. The phase portrait illustrating these equilibrium states 
is plotted in Fig. 1. The separatrix AB, connecting these two singular points, corresponds to 
the "overfall"- type solution with the amplitude E e . With the boundary conditions E = E e 
at z — > — oo, this equation can be written as the equation for wave velocity v e : 



1 - QqV e = X (v E ) = Qev 2 e , 
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and Hq is a constant parameter. Common points of the Xq(ve) function and the quadratic 
parabola Q V E 2 yield the velocity v E we are interested in: 
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The corresponding plot is depicted in Fig. 2. Only the positive value v + is physically 
meaningful. In the limit H — > oo (Bean's model [4]), the wave velocity is defined by 
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In the reverse case, when Hq 
of the wave amplitude E e : 

ve ? 



0, the velocity is largely small and is almost independent 
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For the field H approaching H C2 the following approximation is valid: 
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In conclusion, our results suggest that taking into account an arbitrary dependence of j c 
n H has a minor effect on the wave profile, since the phase portrait remains nearly the 
same. The qualitative analysis presented above may be useful for experimentally studying 
the penetration of a magnetic flux into a superconductor placed into an external magnetic 
field. Improving the protection of superconducting devices at the design stage is another 
possible application of our results. 
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